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Abstract 

We present an exponential F(R) modified gravity model in the Jordan and the 
Einstein frame. We use a general approach in order to investigate and demonstrate 
the viability of the model. Apart from the general features that this models has, which 
actually render it viable at a first step, we address the issues of finite time singularities, 
Newton's law corrections and the scalaron mass. As we will evince, the model passes 
these latter two tests successfully and also has no finite time singularities, a feature 
inherent to other well studied exponential models. 



iff : Introduction 

General Relativity is one of the most sound scientific descriptions of nature, regarding 
gravitational interactions in scales where gravity becomes important. Hence, it has become 
a powerful tool to describe local gravitational interactions and moreover to describe how 
the universe evolves as a whole. The astrophysical observations of the late 90's have put 
into a different perspective the way that we think the universe is evolving. According to 
these observations, the universe has undergone two accelerating phases. The first was the 
inflation period, while the second is the present epoch's acceleration. In reference to the 
latter, this cosmic acceleration is know as dark energy. The observational data for the 
present epoch (the new results of Planck telescope) suggest that the universe is described 
by the ACDM model. The main features of this model is that the universe is almost 
spatially flat, and consists of ordinary matter (~ 4.9%), cold dark matter (~ 26.8%) 
and dark energy (~ 68.3%), with the latter being the reason behind the present day 
acceleration. 

One of the most recent attempts to successfully describe the dark energy is provided by 
the F(R) modified theories of gravity (for informative reviews and very important papers 
on these theories see [IHI]), in the context of which, what changes drastically is not the 
left hand side of the Einstein equations, but the right hand side. In order the Friedmann- 
Robertson- Walker equations give an accelerating solution, the energy momentum tensor 
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must contain a fluid with negative w, and this is achieved with these F(R) theories. Par- 
ticularly, early time inflation and the late time acceleration enjoy a unified description 
within the self-consistent theoretical framework some of these models provide [IHZ1- I n 
addition, some of these models are put into a string theory oriented context [5] and also 
quantum corrections have been calculated [9]. Connections to other gravitational exten- 
sions of general relativity, such as Gauss-Bonnet gravity, have also been pointed out in 
the literature [10] and moreover some exact solutions in strong gravitational background 
have been done in |llj . 

A very interesting and novel feature that most of the F(R) theories of gravity have, is the 
appearance of finite time singularities in the physical parameters. There are several types 
of singularities depending on which physical parameters become singular. This finite type 
singularities can completely change our universe's evolutionary history, leading to rather 
dramatic final states, known as Big Rip, or Little Rip etc. For an important stream of 
papers in regard to cosmological singularities and related issues, see [HE] arid references 
therein. 

With the F(R) theories being a a generalization of general relativity, inevitably, these 
theories are confronted with the successes of general relativity. This entails a number of 
constraints that need to be satisfied, in order a promising modified gravity theory can be 
considered viable. These constraints are related to the local tests of general relativity and 
additionally to various cosmological bounds. The local tests are related to planetary and 
star formation tests. Moreover, each F(R) theory formulated in the Jordan frame has a 
corresponding scalar-tensor gravitational theory in the Einstein frame, the scalarons of 
which have to be classical, so that the theory is quantum-mechanically stable. 

In this present letter, we shall analyze in detail one exponential model which has some 
very appealing features in reference to the constraints we mentioned above. Exponential 
models have been thoroughly studied in references 0IHI7]. One °^ the mos t appealing 
features of the exponential models is that, in some of them, singularities are absent, while 
at the same time these models successfully pass all the tests these models are confronted 
with. The model we shall present has very interesting attributes, since it is free of singu- 
larities. Also we address the matter instability and Newton's law corrections issues and 
as we will demonstrate, the model at hand passes these tests successfully. 
This letter is organized as follows: In section 1 we briefly present the general features of 
F(R) theories in the Jordan frame within the context of the metric formalism. In section 2 
we study in detail the exponential model and also we investigate all the criteria which have 
to be satisfied in order a modified gravity model can be considered viable. In section 3 we 
present briefly an interesting functional resemblance of two F(R) models (one of which is 
the one presented in this article), to the fermi distributions connected to Woods-Saxons 
potentials. The conclusions follow in the end of the article. 



1 General Features of F(R) Dark Energy Models in the Jor- 
dan Frame 

In this section we shall give a brief description of the main features of F(R) theories in 
the Jordan frame. For a much more detailed account consult references [HE]. The F(R) 
modified gravity theories are described by the following four dimensional action: 

S = ^jd 4 xV=gF(R) + S m ( gflI/ ,y m ), (1) 

with k? = 8ttG and S m the matter action of the matter fields "if m . We shall focus on the 
metric formalism for our study and in addition we shall assume that the form of the F(R) 
theory that we shall present is of the form F(R) = R + f(R)- Within the context of the 
metric formalism, by varying the action ([TJ with respect to gw, we obtain the following 
equations of motion: 

F'(R)R^(g) - ^F(R) 9flv - V^J V F'{R) + g^UF'(R) = k% u . (2) 

In the above equation, F'(R) = dF(R)/dR and T^ v is the energy momentum tensor. The 
main idea behind the F(R) modified gravity theories is that, what is actually modified is 
not the left hand side of the Einstein Equations, but the right. Indeed, the above equations 
of motion can be cast in the following form: 



R„u 2 Rg, " u ~ F'(R)V^ + k 



F(R) - RF'(R) 



g»u + V^VuF'(R) - g^UF'{R) 



(3) 
Thus, the energy momentum tensor has another contribution coming from the term: 
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F{R) - RF'(R) 



g^u + V^F'tfl) - g^nF'(R) , (4) 



and this term is what actually models the dark energy. Taking the trace of equation @ 
we obtain the following equation: 

3DF'(R) + RF'(R) - 2F(R) = k 2 T, (5) 

with T being the trace of the energy momentum tensor T = g^T^ u = —p+P, and p, P are 
the energy density and pressure of the matter. This equation actually shows us that there 
is another degree of freedom underlying the F(R) theories, materialized by the scalar field 
F'(R). Consequently, equation ([S} is actually the equation of motion of this scalar degree 
of freedom which is called the "scalaron" . In a flat Friedmann-Lemaitre- Robertson- Walker 
spacetime, the Ricci scalar is given by: 

R = 6(2H 2 + H), (6) 



with H the Hubble parameter, and the "dot" represents differentiation with respect to 
time. Accordingly, the cosmological dynamics are governed by the following equations: 

3F>(R)H* = K\p m + Pr ) + WWR- F(R)) _ 3H p, {Rl (?) 

- 2F'{R)H = K 2 {p m + 4/3p r ) + FF'(R) - HF'(R), 

where p r and p m stand for the radiation and matter energy energy density respectively. 
Hence, the total effective energy density and pressure of matter and geometry are [US]: 
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p + ^[F'{R)R-F{R)-6HF'(R)) (8) 



p + \ ( - F'(R)R + F(R) + 4HF'(R) + 2F'(R) 



F'{R) 

with p, P the total matter energy density and pressure respectively. 

2 General Study and Detailed Study of the Exponential 
Model 

In principle, every viable F(R) model in the metric formalism has to satisfy certain con- 
ditions, in order to be compatible with observations and also consistent with theoretical 
predictions. In this section we shall present a viable exponential model, which has as lim- 
iting cases other exponential models that exist in the literature, and we shall thoroughly 
study the quantitative features of the F(R) theory it describes. The exponential model is 
of the following form, 

F{R) = R - A + B e -^ + ATB' (9) 

with A, B, C, D constants. In the rest of this paper, we shall denote f(R), the second term 
on the right hand side of Eq. ([9]), that is: 

f^-- A + Be-^ + ATB- (10) 

2.1 General Features of the Model 

Let the Ricci scalar of the universe at the present epoch be denoted as Rq ~ 10 -66 eV and 
additionally denote as A, Aj the cosmological constant of the universe today and during 
the inflation period of the universe, respectively. The general features of a viable F(R) 
model are the following [I]: 

• (i) F'(R) > for R > Rq. Moreover if R\ is a de Sitter final attractor of the system, 
then this inequality has to hold true for R> R\. 

• (ii) F"{R) > for R > Rq. This restriction is required for consistency with local 
gravity tests, for the presence of a matter dominated epoch and for the stability of 
cosmological perturbations. 



• (Hi) F(R) = R — Aj for R — > oo. This is required in order inflation occurs in the 
universe. 

• (iv) F(R) = R — A for R—^Rq and f(Ro) = — A. This is required in order the late 
time acceleration occurs. 

• (y) -F(O) = 0, in order a flat spacetime solution exists. 

In order model © satisfies the above constraints, the constants have to satisfy some 
conditions. In order to have a clear picture of how the derivatives behave as a function of 
R, we quote them below: 



F"(R) 



BCe% ( - B + Ae% 



D 2 [B + AeD 
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Due to the square in the denominator of both the derivatives F'(R) and F"(R), it is not 
very difficult to choose the constants A, B, C, D. Indeed, by looking both expressions, 
the first and second derivative are always positive when A > B and D > C. In order 
not to fall into inconsistencies, we also require that the inequality means that the related 
constants have at least ~ 10% difference in their scale, which for example means that A 
is ten times larger than B. Further restrictions on the constants shall be imposed later 
on in this section. Hence, the conditions (i) and (ii) in the list above are satisfied for all 
the values of the Ricci scalar R, an argument that is further supported from the fact that 
the denominator behaves as ~ e 2R ' D , rendering the fraction smaller than one for all R, 
in reference to the first derivative. The second derivative is always positive due to the 
fact that the term — B + Aeo is always positive for A > B due to the exponential term. 
Condition (z/) is automatically satisfied for the model ([9]). 
When R — > oo, the function F(R) behaves as: 

F(R)~R- CB (12) 

A(A + B\ 



It is obvious that A/, the early time cosmological constant that governs the inflationary 
expansion of the universe, must be Aj ~ a(A+b) • Consequently, since Aj ~ 10 20_38 eV 2 , 
the values of A, B, C have to be chosen appropriately in order this requirement is satisfied. 
A set of values that will prove to be appropriate, in order the present and other constraints 
that will be imposed later on in this section are satisfied, are the following: 

A = 0.1, 5 = 0.01, D = 10 22 , C = 10 21 . (13) 

Using the values (|13|) . the inflation period cosmological constant is approximately Aj ~ 
9.09 x 10 20 . After checking these very general features, we proceed to more elaborate 
criteria for the viability of the present exponential model. 



As a final comment before we proceed, notice that when R ^> Rq, the model ([9]), is 
approximately equal to: 

Therefore we can see that the model of [3113] and the model ([9]) coincide in the large cur- 
vature limit. Therefore we expect that these models may give rise to similar cosmological 
dynamics. We shall elaborate on this issue later on in this section. 

2.2 Scalaron Mass and Effective Potential 

The approach we adopt in the present section is based on references [3ll3| . Following OS] , 
the scalaron dynamical degree of freedom is governed by the trace equation of motion ([5]) . 
By making the substitution F'(R) = 1 + f'(R) = e~ x , and performing a perturbation 
around a constant scalar curvature solution R*, so that R = R* + 5R, the equation of 
motion of the scalaron field read: 

™ x -3(-r(^r-^^ = -6(i + / W ) r ' (15) 

with 5x being related to 5R as follows: 

The constant scalar curvature solution R*, is a solution to the equation: 

R* + 2f(R*)-R*f'(R*) = 0. (17) 

The mass square of the scalaron plays an important role in reference to local and planetary 
test of the modified gravity theory. This is equal to: 

Only a positive and large value of the mass square is required, in order no tachyonic 
instability occurs and the corrections to the Newton's law are small. Substituting the 
model (fTUj) . in Eq. (TTTJ) , we have that R* is a de Sitter solution of the system. The value 
R* = corresponds to the Minkowski spacetime, so we would like to see if Minkowski 
spacetime is stable. Particularly, we shall examine if the scalaron theory perturbed around 
the Minkowski solution R* = is stable, or if it has tachyonic or any other sort of 
instabilities. This test is critical for flat spacetime local tests. For R* = and using the 
values (|13p . the mass of the scalaron field for the model (|10p . is equal to: 



M 2 (0) ~ 4.552 x 10 22 . (19) 

The above value is positive and consequently there is no tachyonic instability. In addition, 
since it is very large, the 8R perturbations at long ranges tend to zero and therefore, 



Newton's law has no considerable corrections. Actually, the mass square is positive for a 
large range of values of the parameters, provided that A > B and D > C. This can be 
easily seen by looking the analytic expression of the mass, for the model at hand: 

(A + B)d( -BC+(A + B) 2 d) 
M2(0 > = i(A-B)BC "• (20) 

Before closing this issue, we shall search for other de Sitter solutions of Eq. (|17|) . We 
are especially interested in a late time de Sitter point. Following the notation of [3J, let 
G(R) = F(R) - RF'(R). Since G"(0) < (see Eq. (HH)), the function becomes negative 
and increases after R = Rq. For R = 0(A), the derivatives of the F(R) function behave 
as: 

F'(R) ~ 1, F"(R) ~ 0. (21) 

The computation of the late time de Sitter point is straightforward, since Eq. (|17p . 
simplifies to the following: 

2F(R) = R. (22) 

Searching a solution around R = A and using the values (|13p for the variables A, B, C, D, 
we easily get numerically, that R = A is actually a late time de Sitter solution. 

2.3 Matter Instability Analysis of the Exponential Model 

We now turn our focus on the matter instability issue, which might occur when the scalar 
curvature is large, in reference to the one corresponding to the present epoch. The scalaron 
equation can be written in the form [3j|I]: 



m y (l + f'(R))R 2(R + f(R)) _ k 2 , 

f»(R) p + 3f"(R) 3f»(R) 6f"(R) 



We perturb the scalar curvature around an Einstein gravity solution R e = — ^- > 0, so 
that R = R e + 5R, and we obtain the equation: 

{-dl + U{R e ))5R + const ~ 0, (24) 

with U(R e ) the effective potential, which is: 



F""(R e ) F'"(R e ) 2 ^ Re 

F"(R e ) F"(R e ) 2 r pReV Re+ 3 



U{Re) = ( ^77^ " W£k VpfleV'*. + 7- (25) 



F'{R e )F"'{R e )R e F'(R e ) + 2F(R e )F'"{R e ) F'"(R e )R e 



3F"{R e ) 2 3F"(R e y 3F"(R e y 3F"(R e ) 2 ' 

In the case the effective potential is positive, the perturbation 5R becomes exponentially 
large and the whole system is rendered unstable [21H]. Therefore, the matter stability 
condition is U(R e ) < 0. Let us see what happens in the case of the exponential model 



at hand ([9|). If we substitute ([9]) to Eq. (|25|) . and keep only the leading order terms, the 
potential U(R) can be written as follows: 



D 2 e-^[B + Ae^ 



U(R e ) j> ^-. (26) 

3BC B - Aeo 



Due to the term in the denominator, namely B — Aeo, and since we imposed the condition 
A > B for the viability of the model, the potential given by relation (j26[) is always negative, 
for all curvature values and for all the values of the parameters so that, A > B and D > C. 
Thereby, the matter stability condition is fulfilled for the model ([9]). 

Before finishing this section we would to address an issue related to the matter era period 
of the cosmological model. Following |4J, during the matter era period, and neglecting the 
contribution of radiation, one has: 

_ p RF'(R) _ 

peff ~r(Ry "W"' ( } 

with p e ff and p the total energy density (matter and dark energy) and p the matter energy 
density. From the second relation of Eq. (|2"Tj) we get [I] : 

tU = °' (28) 

So at the matter era we must have that F"(R) ~ 0. As pointed in [JJ, an .F(i2)-theory is 
acceptable if the modified gravity contribution vanishes during this era and also F'(R) ~ 1. 
In addition, the second derivative must be very small and positive [3J. Let us examine what 
happens for the model ([9]). As can be seen, the second derivative of the F(R) exponential 
is governed by the exponential term in the denominator, that is: 

F"(R) - -££*• (29) 

Using the values (fT3j) for the parameters A, B,C, D and due to the exponential dependence, 
it can be easily checked that the second derivative during the matter era is very close to 
zero and positive. 

2.4 Einstein Frame Analysis of the Model 

Apart from the latter approach we adopted in order to address the problem of matter 
instability, there is another approach for this problem which is related to the Einstein 
frame of the F(R) theory. Hence, we shall investigate the problem of matter instability 
and Newton's law corrections, in the Einstein frame. The Jordan frame and the Einstein 
frame description of F(R) modified theories of gravity are mathematically equivalent. 
Actually, the F(R) theories in the Jordan frame become scalar tensor theories with a 



potential in the Einstein frame [UG2E1E] • Following reference [3], an auxiliary field A is 
introduced in the action ([I]) 

S = ±fdx*y/=j((R-A)(l + f(A))+A + f(AJ). (30) 

Using the conformal transformation g^ — > e a g^u, with a = — ln(l + f'{A)), and bearing 
in mind that the equations of motion with respect to .4, yield the result that A = R, we 
obtain the Einstein frame action: 

Se = ^J dx 4 V=g(R - \g» v d»od v o - V(a)) , (31) 

with the effective potential V{a), being equal to: 



,2«rw-/--^_ A F ^ A ) 



V(a) = e° g (e-°) - e»F{g(e-)) = ^ - -^ (32) 

The mass corresponding to the field a is equal to: 

° 2 da 2 ~ 2\F'(A) F'{A) 2 F"{A))' { } 

In order the Newton's law corrections are small, this mass term has to be a large number 
and of course positive in order tachyonic instabilities are avoided. Let us see how this 
term behaves. Substituting Eq. ([9]) into (j33j) . the mass of the a field is equal to (we keep 
only dominating terms): 

< ^ -^ (34) 

Let us find the value of the above mass, by using the values (|13p for the parameters. For a 
wide range of values of the scalar curvature, the mass (|34p has the value m 2 ~ 7.2 x 10 22 , 
which is large, and thus we conclude that the Newton's law corrections are negligible. In 
addition, since it is positive, the theory is free from tachyonic instabilities. 

2.5 Analysis of Finite Time Singularities of the Exponential Model 

In this section we shall analyze in detail the singularity structure of the model at hand. 
The singularities that quite frequently occur in F(R) theories of gravity have the form 
H{t) = h/(to — t)P , with h and to positive constants [4j. There are four types of finite 
time future singularities [HHIE] which are listed below: 

• Type I (Big Rip): For t — > to, a(t) — } oo, p e jf — > oo and \p e ff\ — ► oo ((3 > 1) 

• Type II (sudden): For t — } to, a(t) — } ao, p e ff — > Po and \p e ff\ —* oo (—1 < /3 < 0) 

• Type III: For t — > to, a(t) — > ao, p e ff — > oo and \p e ff\ — S- oo (0 < /3 < 1) 

• Type IV: For t — > to, a(t) — > ao, p e ff — > and \p e ff\ —* and higher derivatives of 
the Hubble parameter H diverge. 



The singularity structure of the F(R) model © is quite similar to the one presented 
in reference [3J. Indeed, at the limit R —> oo, we have: 

lim F(R) ~ R + const, lim F'(R) ~ 1. (35) 

In addition, the higher order derivatives of the F(R) function fl5J), tend to zero exponen- 
tially. Therefore, neither Type I nor Type III singularities can appear in the present model, 
which is exactly what happens in reference [5J . Now let us consider the other two possible 
singularities, namely Type II and Type IV. In reference to the first, when R — > — oo, the 
F(R) function can be written as: 

F{R)^R-^e- R ' D + ^. (36) 

Consequently, the total effective energy density p e ff of Eq. (JSJ), exponentially decays 
for the model ([!]) and also the total effective pressure behaves analogously. Thereby, the 
Type II singularity cannot occur, since the pressure does not diverge. In addition, Type IV 
singularities cannot be realized. The reason is in absolute concordance with the arguments 
of [3J. Indeed, when R—tO, the F(R) function behaves as 

F(R ^ R -<ATWd R - (37) 

and therefore, the effective energy density behaves as ~ l/(*o — t) an d is larger than 
1/(^0 — t) 2 ^ j f° r ft < — 1- Hence the Type IV singularity cannot occur. As in the [3J case, 
the model we presented is free of singularities. The only difference between the two models 
is the different reasoning for the Type II singularity, since in the present model it is the 
finite pressure that renders the model free of these singularities, in contrast to the model 
of [JJ, where the effective energy density diverges. In addition, since according to Eq. 
(|14|). the two models overlap when R — > oo, the reasoning on why the Type I and Type II 
singularities do not occur in both models is identical in both cases, as it was expected. 

3 Functional Resemblance of two F(R) Models to Wood- 
Saxons Nuclear Potentials 

It worths to note that the model we used in this letter described by the F(R) function 
§§§, has many similarities to some Woods-Saxons potentials that are used to describe the 
Nucleon interactions. The Woods-Saxon potentials are phenomenological potential for 
the nucleons in the atomic nucleus. It is used in the shell model and describes the forces 
applied to each nucleon. The general form of the potential is, 

V(r) = \,, (38) 

1 + e a 

where r the distance from the center of the nucleus. Note one particular fermi function 
distribution, related to (f38|) 

u / (r)~ J B c tanh(^). (39) 
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Note that this function ([39]) has great similarities to the Tsujikawa model [IJ[2], for F(R) 
modified gravitational theories, namely: 

F{R) = R-R C tanh(R/R c ) . (40) 

Indeed, if we change r — > —R, and add R in expression (f3T^ we get the Tsujikawa model. 
Using the same line of reasoning, the following Fermi shape distribution, 

"/CO = T— ~ r , (41) 

leads to the model we used in this article. Of course there is no physical connection 
between the two systems, the nucleus potential and the gravitational systems. We wanted 
to stress the fact that by performing the same transformations to two different nuclear 
fermi distributions, we get two viable F(R) dark energy models. 

Concluding Remarks 

In this article, we studied in detail an exponential model of F(R) gravity in the metric 
formalism. Apart from the very general features that make the model viable at a first 
step, the model has also other interesting features that render it promising. Particularly, 
the scalaron mass is positive and large, which means that the theory is free of tachyonic 
instabilities and free of matter instabilities. The matter instability and Newton's law 
corrections problem has been addressed in both the Jordan and Einstein frames, and 
the result verified once more that the model does not lead to problems or inconsistencies. 
Moreover, we studied the finite time singularities issue and as we demonstrated, the model 
is free of the four types of finite time singularities. In addition, we also presented what 
motivated us to use such an exponential model. As we showed, this model and another 
F(R) model have functional resemblance to two Fermi distributions corresponding to 
Woods-Saxons potentials. 

We intended to show the basic features of this F(R) model, using a very general approach, 
similar to the one used in the articles [3JII]. It is interesting to address other interesting 
issues, such as oscillations of the dark energy, or curvature singularities in dense gravita- 
tional backgrounds |12j . which however are beyond the scope of this letter. We hope to 
address such issues in a future work. 

References 

[1] Shinji Tsujikawa, Lect. Notes Phys. 800 (2010) 99; Shin'ichi Nojiri, Sergei D. 
Odintsov, Int. J. Geom. Meth. Mod.Phys. 4 (2007) 115; M.C.B. Abdalla, Shin'ichi 
Nojiri, Sergei D. Odintsov, Class. Quant. Grav. 22 (2005) L35; Antonio De Felice, 
Shinji Tsujikawa, Living Rev. Rel. 13 (2010) 3; Thomas P. Sotiriou, Class. Quant. Grav. 
26 (2009) 152001; Thomas P. Sotiriou, Valerio Faraoni, Rev.Mod.Phys. 82 (2010) 
451-497; Shin'ichi Nojiri, Sergei D. Odintsov, Phys.Rept. 505 (2011) 59; M. Akbar, 



11 



Rong-Gen Cai, Phys.Lett. B635 (2006) 7-10; Marek Szydlowski, Aleksandra Kurek, 
Int.J.Geom.Meth.Mod.Phys. 4 (2007) 3; Andrzej Borowiec, Wlodzimierz Godlowski, 
Marek Szydlowski, Int.J.Geom.Meth.Mod.Phys. 4 (2007) 183 

[2] Shin'ichi Nojiri, Sergei D. Odintsov, Phys.Rev. D70 (2004) 103522; Thomas P. 
Sotiriou, Class. Quant. Grav. 23 (2006) 5117; Thomas Faulkner, Max Tegmark, Emory 
F. Bunn, Yi Mao, Phys.Rev. D76 (2007) 063505; Shinji Tsujikawa, Kotub Uddin, 
Shuntaro Mizuno, Reza Tavakol, Jun'ichi Yokoyama, Phys.Rev. D77 (2008) 103009; 
Shinji Tsujikawa, Phys.Rev. D77 (2008) 023507; Shinji Tsujikawa, Phys.Rev. D76 
(2007) 023514; Shin'ichi Nojiri, Sergei D. Odintsov, J.Phys. A40 (2007) 6725-6732; 
Shin'ichi Nojiri, Sergei D. Odintsov, Hrvoje Stefancic, Phys.Rev. D74 (2006) 086009; 
Shin'ichi Nojiri, Sergei D. Odintsov, Phys.Rev. D74 (2006) 086005; Shin'ichi No- 
jiri, Sergei D. Odintsov, Phys.Lett. B595 (2004) 1; Emilio Elizalde, Shin'ichi Nojiri, 
Sergei D. Odintsov, Phys.Rev. D70 (2004) 043539; Luca Amendola, Shinji Tsujikawa, 
Phys.Lett. B660 (2008) 125; Stephane Fay, Reza Tavakol, Shinji Tsujikawa, Phys.Rev. 
D75 (2007) 063509; Luca Amendola, Radouane Gannouji, David Polarski, Shinji 
Tsujikawa; Phys.Rev. D75 (2007) 083504; Luca Amendola, David Polarski, Shinji 
Tsujikawa, Phys.Rev.Lett. 98 (2007) 131302; Shin'ichi Nojiri, Sergei D. Odintsov, 
Shinji Tsujikawa Phys.Rev. D71 (2005) 063004; Kazuharu Bamba, Shin'ichi Nojiri, 
Sergei D. Odintsov, JCAP 0810 (2008) 045; Shin'ichi Nojiri, Sergei D. Odintsov 
Phys.Rev. D78 (2008) 046006; Shin'ichi Nojiri, Sergei D. Odintsov, Phys.Lett. B599 
(2004) 137; S. Capozziello, S. Nojiri, S.D. Odintsov, Phys.Lett. B632 (2006) 597; 
Shin'ichi Nojiri, Sergei D. Odintsov Gen.Rel.Grav. 38 (2006) 1285-1304; Shin'ichi No- 
jiri, Sergei D. Odintsov, Diego Saez-GomezAug 2009. 11 pp. Phys.Lett. B681 (2009) 
74; Kazuharu Bamba, Antonio Lopez-Revelles, R. Myrzakulov, S.D. Odintsov, L. 
Sebastiani, Class.Quant.Grav. 30 (2013) 015008; E. Elizalde, A.N. Makarenko, S. 
Nojiri, V.V. Obukhov, S.D. Odintsov, Astrophys. Space Sci. 344 (2013) 479-488; Pe- 
ter K.S. Dunsby, Emilo Elizalde, Rituparno Goswami, Sergei Odintsov, Diego Saez 
Gomez, Phys.Rev. D82 (2010) 023519; Masud Chaichian, Shin'ichi Nojiri, Sergei D. 
Odintsov, Markku Oksanen , Anca Tureanu Class.Quant.Grav. 27 (2010) 185021 

[3] G. Cognola, E. Elizalde, S. Nojiri, S.D. Odintsov, L. Sebastiani, S. Zerbini Phys.Rev. 
D77 (2008) 046009 

[4] E. Elizalde, S. Nojiri, S.D. Odintsov, L. Sebastiani, S. Zerbini, Phys.Rev. D83 (2011) 
086006 

[5] Shin'ichi Nojiri, Sergei D. Odintsov Phys.Lett. B716 (2012) 377; Kazuharu Bamba, 
Shinichi Nojiri, Sergei D. Odintsov Gen.Rel.Grav. 44 (2012) 1321; G. Cog- 
nola, E. Elizalde, S.D. Odintsov, P. Tretyakov, S. Zerbini Phys.Rev. D79 (2009) 
044001;Kazuharu Bamba, Chao-Qiang Geng, Shin'ichi Nojiri, Sergei D. Odintsov 
Phys.Rev. D79 (2009) 083014; S. Jhingan, S. Nojiri, S.D. Odintsov, M. Sami, I 
Thongkool, S. Zerbini Phys.Lett. B663 (2008) 424; Emilio Elizalde, Shin'ichi No- 
jiri, Sergei D. Odintsov, Diego Saez-Gomez, Valerio Faraoni Phys.Rev. D77 (2008) 
106005; Shin'ichi Nojiri, Sergei D. Odintsov Phys.Lett. B659 (2008) 821; Shin'ichi 



12 



Nojiri (Nagoya U.), Sergei D. Odintsov Phys.Lett. B657 (2007) 238; Shin'ichi No- 
jiri, Sergei D. Odintsov Phys.Lett. B639 (2006) 144; Salvatore Capozziello, S. No- 
jiri, S.D. Odintsov, A. TroisiPhys.Lett. B639 (2006) 135; Kazuharu Bamba, Chao- 
Qiang Geng, Shin'ichi Nojiri, Sergei D. Odintsov, Mod.Phys.Lett. A25 (2010) 900; 
Artyom V. Astashenok, Sergei D. Odintsov Phys.Lett. B718 (2013) 1194; Sergio 
E. Joras, Int.J.Mod.Phys. A26 (2011) 3730; Yuki Yokokura, Int.J.Mod.Phys. A27 
(2012) 1250160; Nikodem J. Poplawski, Int.J.Mod.Phys. A23 (2008) 1891; Salvatore 
Capozziello, Christian Corda, Mariafelicia De Laurentis Mod.Phys.Lett. A22 (2007) 
1097 

[6] F. Briscese, E. Elizalde, S. Nojiri, S.D. Odintsov, Phys.Lett. B646 (2007) 105; E. 
Barausse, T.P. Sotiriou, J.C. Miller, Class.Quant.Grav. 25 (2008) 105008; Shin'ichi 
Nojiri, Sergei D. Odintsov Gen.Rel.Grav. 36 (2004) 1765; Shin'ichi Nojiri, Sergei 
D. Odintsov Phys.Rev. D68 (2003) 123512; Shin'ichi Nojiri, Sergei D. Odintsov 
Phys.Rev. D72 (2005) 023003; Gianluca Allemandi, Andrzej Borowiec, Mauro Fran- 
caviglia, Sergei D. Odintsov Phys.Rev. D72 (2005) 063505; Artyom V. Astashenok, 
Sergei D. Odintsov Phys.Lett. B718 (2013) 1194; Artyom V. Astashenok, Shinichi No- 
jiri, Sergei D. Odintsov Phys.Lett. B709 (2012) 396; Paul H. Frampton, Kevin J. Lud- 
wick, Shin'ichi Nojiri, Sergei D. Odintsov, Robert J. Scherrer Phys.Lett. B708 (2012) 
204; I. Brevik, E. Elizalde, S. Nojiri, S.D. Odintsov Phys.Rev. D84 (2011) 103508; 
Shin'ichi Nojiri, Sergei D. Odintsov, Alexey Toporensky Gen.Rel.Grav. 42 (2010) 
1997; Shin'ichi Nojiri, Sergei D. Odintsov Phys.Lett. B686 (2010) 44; S. Capozziello, 
M. De Laurentis, S. Nojiri, S.D. Odintsov Phys.Rev. D79 (2009) 124007; Kazuharu 
Bamba , Sergei D. Odintsov Phys.Rev. D77 (2008) 123532 

[7] K. Bamba, Chao-Qiang Geng, Chung-Chi Lee, JCAP, 08 (2010) 021; E. Elizalde, S.D. 
Odintsov, L. Sebastiani, S. Zerbini, Eur.Phys.J. C72 (2012) 1843; Eric V. Linder, 
Phys.Rev. D80 (2009) 123528; Louis Yang, Chung-Chi Lee, Ling- Wei Luo, Chao- 
Qiang Geng, Phys.Rev. D82 (2010) 103515; Pengjie Zhang; Phys.Rev. D73 (2006) 
123504 

[8] Emilio Elizalde, James E. Lidsey, Shinichi Nojiri, Sergei D. Odintsov, Phys.Lett. B574 
(2003) 1; Shin'ichi Nojiri, Sergei D. Odintsov Phys.Lett. B576 (2003) 5; Shin'ichi 
Nojiri, Sergei D. Odintsov, M. Sami Phys.Rev. D74 (2006) 046004; E. Elizalde, S. 
Jhingan, S. Nojiri, S.D. Odintsov, M. Sami, I. Thongkool, Eur.Phys.J. C53 (2008) 
447; Guido Cognola, Emilio Elizalde, Shin'ichi Nojiri, Sergei Odintsov, Sergio Zerbini, 
Phys. Rev. D75 (2007) 086002 

[9] Shin'ichi Nojiri, Sergei D. Odintsov, Phys.Lett. B562 (2003) 147; Shin'ichi No- 
jiri, Sergei D. Odintsov, Mod.Phys.Lett. A19 (2004) 1273; Shin'ichi Nojiri, Sergei 
D. Odintsov, Mod.Phys.Lett. A19 (2004) 627; Shin'ichi Nojiri, Sergei D. Odintsov 
Phys.Lett. B565 (2003) 1; Shin'ichi Nojiri, Sergei D. Odintsov, Phys.Lett. B637 (2006) 
139; Kazuharu Bamba, Jaume de Haro, Sergei D. Odintsov, JCAP 1302 (2013) 008; 
Kazuharu Bamba, Shin'ichi Nojiri, Sergei D. Odintsov, Phys.Rev. D85 (2012) 044012; 
Emilio Elizalde, Shin'ichi Nojiri, Sergei D. Odintsov, Peng Wang, Phys.Rev. D71 



13 



(2005) 103504; Guido Cognola, Emilio Elizalde, Shin'ichi Nojiri, Sergei D. Odintsov, 
Sergio Zerbini, JCAP 0502 (2005) 010; Guido Cognola, Sergio Zerbini, J.Phys. A39 

(2006) 6245; T. Tahamtan, O. Gurtug, Eur.Phys.J. C72 (2012) 2091; Ali Shojai, 
Fatimah Shojai, Gen.Rel.Grav. 40 (2008) 1967 

[10] Shin'ichi Nojiri, Sergei D. Odintsov, Phys.Lett. B631 (2005) 1; Guido Cognola, Emilio 
Elizalde, Shin'ichi Nojiri, Sergei D. Odintsov, Sergio Zerbini, Phys.Rev. D73 (2006) 
084007; Shin'ichi Nojiri, Sergei D. Odintsov, Misao Sasaki, Phys.Rev. D71 (2005) 
123509; Shin'ichi Nojiri, Sergei D. Odintsov, Petr V. Tretyakov, Phys.Lett. B651 

(2007) 224 

[11] Alireza Hojjati, Levon Pogosian, Alessandra Silvestri, Starla Talbot, Phys.Rev. D86 
(2012) 123503; J.P. Morais Graca, V.B. Bezerra, Mod.Phys.Lett. A27 (2012) 1250178; 
M. Sharif, Sadia Arif, Mod.Phys.Lett. A27 (2012) 1250138; H. Saiedi, Mod.Phys.Lett. 
A27 (2012) 1250220; T.R.P. Carames, E.R. Bezerra de Mello, M.E.X. Guimaraes, 
Mod.Phys.Lett. A27 (2012) 1250177; Solmaz Asgari, Reza Saffari, Gen.Rel.Grav. 44 
(2012) 737; K.A. Bronnikov, M.V. Skvortsova, A. A. Starobinsky, Grav.Cosmol. 16 
(2010) 216; A. Azadi, D. Momeni, M. Nouri-Zonoz, Phys.Lett. B670 (2008) 210; 
Reza Saffari, Sohrab Rahvar, Mod.Phys.Lett. A24 (2009) 305-309 

[12] E.V. Arbuzova, A.D. Dolgov, Phys. Lett. B 700 (2011) 289; Kazuharu Bamba, 
Shin'ichi Nojiri, Sergei D. Odintsov, Phys.Lett. B 698 (2011) 451; Chung-Chi Lee, 
Chao-Qiang Geng, Louis Yang, Prog. Theor. Phys. 128 (2012) 415 



14 



